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Hartree-Fock model is presented. The objective is to demonstrate that the terminating states, 
due to their intrinsic simplicity, offer unique and so far unexplored opportunities to study different 
aspects of the effective NN interaction or nuclear local energy density functional. In particular, 



^ \ we demonstrate that the agreement of the calculations to the data depend on the spin fields and 
00 ' 

' the spin-orbit term which, in turn, allows to constrain the appropriate Landau parameters and the 

CD ' strength of the spin-orbit potential. 
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I. INTRODUCTION 



Within the mean-field approach the spin-orbit {is) sphtting is usually studied via the 
comparison of theoretical and experimental single-particle {sp) energies of the £s-doublet. 
The requirement is that both £s-partners should be simultaneously occupied [l|, 0|. The 
method assumes that under such conditions the core polarization effects, which are known 
to modify strongly sp energies jsl 0] , are similar for £s-partners and therefore do not affect 
the £s-splittings, at least not in a major way. 

The method requires by definition precise empirical knowledge on sp energies of deep-hole 
states which are very difficult to measure. In addition, particle vibration coupling may con- 
tribute to the splitting and perturb the pure sp picture. The available data on ^s-splittings 
and their isotopic or isotonic dependence are therefore both scarce and uncertain. For ex- 
ample, in the A~40 mass region, which is of primary interest in this work, the most recent 
evaluations give Aed^^^^d5/2^^ MeV in "^^Ca and Ae^g^^-rfs/a^ 5 MeV in "^^Ca, respectively, 
while older works give Aedy^-d./^^Q-SMeV 0, 7.3 MeV Q, and ^7.7 MeV |^ in ^°Ca 
and ^5.3 MeV [s^ in ^^Ca. More detailed information on sp levels can be found in Ref. j^. 

In this work we want to pursue a novel method to study the £s-potential by exploring 
high-spin data. The method is based on a direct comparison of the excitation energies 
of terminating states [which are maximum-spin states within a given sp configuration] for 
two carefully selected configurations. In this study we limit ourselves to the '^3/2/7/2^ 
fy2 configurations in A~50, 20^Z<N^24 nuclei, where n denotes the number of valence 
particles outside the ^''Ca core. The difference, AE, between the excitation energies of 
states terminating within the 1^3/2/7/2^ fy^ configurations is dominated by the size 
of the magic gap 20. The magnitude of the magic gap 20 is in turn governed by the 
strength of the £s-potential. Indeed, for the spherical Nilsson Hamiltonian i.e. the 
three-dimensional harmonic-oscillator (HO) potential augmented by a one-body spin-orbit 
—2nfujJols and orbit-orbit —Kfihuoi'^ term one has: 
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Hmisson - -hujo = hiUo{N-K [2£s + _ {£'^)n)] } . (1) 

Hence, within the Nilsson model, which is considered to be a fundamental approximation 
for the nuclear mean-field potential, the magnitude of the magic gap 20, or more precisely 
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the /7/2— c?3/2 splitting reads: 

Ae2o = ^0(1 — 6k — 2K/i). (2) 

In light nuclei the nuclear potential resembles the pure HO, thus /i ~ 0. Hence, within the 
Nilsson model, the is-teim plays a dominant role in establishing the size of the magic gap 20. 
Indeed, hwo determines the global energy scale in low energy nuclear physics and one does 
not expect this value to change substantially. Although Eq. pertains to the /7/2— '^3/2 
splitting, the conclusion drawn above seems to be easily extendable to heavier nuclei since 
there /i^l/2 as a consequence of approximate pseudo-SU(3) symmetry [ll, 12, 13 1. 

By limiting ourselves to the study of band terminating states only, we access the regime of 
essentially unperturbed sp motion, where correlations going beyond mean-field are expected 
to be strongly suppressed. Indeed, the success of simple Nilsson-Strutinsky calculations 
of terminating bands by Ragnarsson and coworkers, for review see Q], nicely confirm the 
structural purity and sp nature of the terminating states. However, within the self-consistent 
approaches, in particular within the Skyrme-Hartree-Fock (shf) approach which is used here, 
there are additional difficulties and in turn uncertainties related to the limited knowledge 
of the time-odd components of the mean-field. Direct studies of these terms are not only 
scarce but also in many cases inconclusive, see and refs. quoted therein. In contrast, 
the structural simplicity of terminating states is appealing and allows for a direct study of 
these terms. 

The method proposed here to determine the is-term has clear advantages over the stan- 
dard method mentioned at the beginning: («) it uses terminating states which are probably 
the best examples of unperturbed single-particle motion; {ii) the terminating states are 
uniquely defined implying that configuration mixing going beyond mean-field is expected 
to be marginal; [iii) shape polarization is included automatically within the calculations 
scheme and no further ad hoc assumptions are necessary; (iv) a wide set of rather precise 
experimental data is already available throughout the periodic table. 

The paper is organized as follows. Sect. ITTl overviews available empirical data. Sect. IIIII 
briefiy recalls the local energy density functional (ledf) formulation of the SHF method. 
Sect. II VI reveals the problems related to spin fields emerging from the Skyrme force induced 
local energy density functional (s-ledf) in N^^Z nuclei both in the ground state as well as 
at the band termination. It is shown, that an unification of spin fields cures these problems 
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Ref. 


•^7/9 ■ -^[^maa:^] ^max 
' / 


J-l fn+l . T^\T 1 T 

'M')J 7 1') ' [-^ max i ^max 




iCa22 [18j 


3.189 6+ 


8.297 11" 


5.108 


!^Ca24 [19j 


10.568 8+ 


5.088 13- 


5.480 


4|SC23 [20] 


9.141 11+ 


3.567 15- 


5.574 


!fSc24 [21J 


5.417 23/2" 


11.022 31/2+ 


5.605 






15.701 35/2~ 


10.284 


iTi23 [19] 


7.143 27/2- 


13.028 33/2+ 


5.885 


iTi24 [1^ 


10.034 14+ 


15.549 17- 


5.515 


23V24 [22J 


10.004 31/2" 


15.259 35/2+ 


5.255 



TABLE I: Spins and excitation energies of terminating states in 20 ^ Z < ^ 24 nuclei. The 
first two columns are representative the for /"^g configuration where n denotes the number of 
valence particles outside the ^'^Ca core. The next two columns are representative for the A'^J'^^^^ 
configuration involving the Ip-lh proton excitation across the magic gap 20. For the case of ^^Sc, 
data on the ^3/2/7/2^ configuration involving two Ip-lh excitations from the to fY/2 are also 
included. The relative excitation energies between the two configurations are given in the last 
column. 



leading to an unified description of terminating states. The remaining discrepancy between 
experiment and theory can further be reduced by tuning (weakening) the strength of the spin- 
orbit interaction as shown in Sect. |3 Finally, in Sect. I VII we test the additivity relation for 
the only known case of a 2p-2h excitation across the magic gap in ^^Sc. All SHF calculations 
presented in this paper were done using the SHF code hfodd of Dobaczewski, Dudek, and 



Olbratowski 



16 



17|. 



II. EXPERIMENTAL DATA 

All available experimental data on the terminating states for fy^ and c?;^/2/7/2^ configu- 
rations in A~50, 20^Z<N^24 nuclei where both states are known are listed in Tab. HJ In 
the present data set we have excluded the N=Z nuclei since in these nuclei the terminating 
state dzLf 7/2^ uniquely defined. Indeed, two states with the same value of Imax can 
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be formed, having approximately (due to the weak isospin symmetry breaking) symmetric 
(T = 0) and antisymmetric (T = 1) proton- neutron configurations, respectively. N=Z 
nuclei will be addressed in a forthcoming paper 

The differences in excitation energies, AE^xp, listed in the last column of Tab. |1] are 
fairly constant. The mean value equals AEexp = 5.489 MeV while the standard deviation 
is 0" = 0.251 MeV i.e. at the level of ~5% only. This suggests that the bulk part of AEexp 
is indeed related to the energy of Ip-lh excitation across the gap 20, and that polarization 
effects are either weak or, most likely, cancel out. 

This conclusion is further supported by a recent measurement in ^^Sc of a terminating 
state at spin J^^^ = 35/2^ and excitation energy of 15.701 MeV This state involves the 
2p-2h excitation from the 71^3/2 to vr/7/2 sub-shell. The excitation energy with respect to the 
aligned fy^ state equals AE2p2h = 10.284 MeV. For the extreme case of pure additivity this 
excitation can be composed from two, uniquely defined, Ip-lh states terminating at I^^^ = 
29/2^ (termination at unfavored signature) which both are known. The empirical excitation 
energies of these states, relative to the aligned state, are 4.753 MeV and 5.783 MeV, 
respectively. The resulting sum yields Ai52p2ft =10.536 MeV implying that additivity holds 
within ~2%. 

III. SKYRME-HARTREE-FOCK LOCAL ENERGY DENSITY FUNCTIONAL 

The starting point of the SHF approach is an energy density functional which, in the 
isoscalar-isovector t = 0,1 representation, takes the following form: 



The local energy density functional Ti. (ledf) is uniquely expressed as a bilinear form of 
time-even (te) p, r, J and time-odd (to) s,T,j local densities, currents, and by their 
derivatives: 




t=o,i 



(3) 



nr\r) = C^pI + Ct'ptApt + Clptn + Ct J I + Cj'ptV ■ J, 



(4) 



Sfiy xjt). 



(5) 
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The division of ledf into te and to parts is very convenient since the latter contributes 
only when time reversal symmetry is broken. 

In the above formulas J f = 'l2f_iu "^fiu^t while the vector spin-orbit density, Jf, is not 
an independent quantity but constitutes an antisymmetric part of the tensor density, i.e., 
•^t = S^!/ ^Ati^'^Ati',*- Definitions of all local densities and currents p, r, J ,s,T,j can be 
found in numerous references and will not be repeated here. We follow the notation used 



in Refs. 



25l | where references to earlier works can be found as well. 



By taking an expectation value of the Skyrme force over the Slater determinant one 
obtains LEDF (jSI)^® with 20 coupling constants C that are expressed uniquely through the 
10 parameters i = 0, 1,2,3 and W, a of the standard Skyrme force. The appropriate 



formulas can be found, e.g., in Refs. 



2J, |25|. Due to t 



le local gauge invariance (which 
26j only 14 coupling constants C 



includes the Galilean invariance) of the Skyrme force 
are independent quantities. The local gauge invariance links three pairs of time-even and 
time-odd constants in the following way: 



(6) 



Since the SHE approach uses interaction-induced coupling constants C it constitutes a 
restricted version of the local energy density theory of Hohenberg-Kohn-Sham j27, 3, Si 



31|, and Sly5 



Other 



type. However, only very few she approaches rigidly enforce the Skyrme- force- related values 
of C. Among those studied here these include SkP js^, SkXc 

forces studied here, including Sly4 32|, SIII 33[, SkO [3J], and SkM* j35| disregard the 
tensor terms by setting C/ = Cj^ = 0. This is not only due to practical reasons (these terms 
are the most difficult technically) but also due to lack of clear experimental information 
that would allow for reasonable estimates of their strengths. Moreover, in the case of SkO 
we were forced to set Cf^^ = to assure convergence. All versions of lede that use the 
Skyrme-force-induced C values, including those taking C/ = = will be called later 
Skyrme-LEDE (s-LEDE). 



IV. THE SPIN FIELDS 



The strengths of the spin fields, C^s'^ and Cf^'^s-As, emerging from the Skyrme force 
appear to be more or less accidental. This is illustrated in Fig. ^ showing the isoscalar 
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FIG. 1: Isoscalar Landau parameters go (solid circles) and gi (open circles) for various Skyrme 
forces. Vertical lines marks the values 50=0. 4 and (7i=-0.19 recommended by Bender et al. 
from their study of the Gamow- Teller resonances. 



Landau parameters qq and gi calculated for the Skyrme forces under study. The Landau 
parameters are related to the ledf strengths in the following way 25^ : 



(7) 



^?o' = Noi2Ct + 2Cff3pT), g[ = -2NoCf/3p'J\ (8) 

where P = (37rV2)2/3, and N^^^ = ■K^h^/2m*kF is an effective-mass-dependent normalization 
factor. The Skyrme-force-induced go and gi parameters, see Fig. d are indeed scattered 
rather randomly reflecting the fact that Skyrme forces are fitted ultimately to the te channel 
while the TO components of the S-LEDF are only cross-checked mostly through the high- 
spin (cranking) applications. In Ref. j3| the preferred values (7o=0.4, (7g=1.2 and gi=- 
.62 have been established from the analysis of the Gamow- Teller resonances, see 



0.19, g[=Q 
also Ref. 



36| and refs. quoted therein. These values were obtained under an additional 
assumption of the density independence of Cf and for C^^=Q. The ledf with spin fields 



defined in this way will be called later the Landau-LEDF (l-ledf). 



A. The spin fields at the ground state 

Before proceeding to the study of the terminating states let us discuss the spin fields of 
the ground states (calculated in HF approximation) of N^Z nuclei. It is known that the 
binding energies calculated using the complete s-ledf functional exhibit a rather peculiar 
behavior in odd-odd Tz = and some odd- A |Tj| = 1/2 nuclei j^^]- It manifests itself via 
additional binding energy of the order of ~ 1 MeV as compared to the SHF calculations using 
only the TE part of the S-LEDF. The effect disappears for ITjI > 1 nuclei. 

The enhancement in the binding energy of N^Z nuclei is due to a strong polarization 
effect exerted by the spin field of the odd-particle (s) on the spin field of the core, as illustrated 
in Figs. 121 and El for a representative example of the manganese isotopes. Fig. |21 shows the 
Oxy component of the spin density s± in ^°Mn for three selected cross sections through 
this axially deformed nucleus, that include the near-equatorial plane at z=0.27fm as well 
as 2;=1.35fm and 2.46 fm planes. To visualize the polarization effect we decompose the spin 
density into contributions of the valence particles (note that ^ in this N=Z nucleus) 
and the core. The topology of the surfaces shown on the left hand side reflects the structure 
of the dominant asymptotic [312]5/2 Nilsson component in the wave function of the valence 
particle. Indeed, the one particle contribution to the spin field in a simplex-conserving axial 
HO basis state: 

^Afn.|A|;s=±i = -^(^7Vn,A;l/2 ± Nn,-A--l/2) whcrC Nn,A = V'A^n. [A| 6*'^'^, (9) 

is 

= -^|^7Vn,|A|(p, 2)Psin2A(y9, Sy = -^\ipNn,\A\{p,z)\^ cos2A(f, Sz = 0, (10) 

see Ref. [l^ for further details. Hence, |sa;|/|sy| = |ctg2A(^| and the plot of s±_ shows the 
characteristic "vortex" lines for = 1, i.e., for ^9 = [45° + r27r]/2A where n = 0,1,2,.... 

Two such lines that appear in Fig. |21are consistent with A = 2. Moreover, the small values 
of Sxy ~ over the entire equatorial plane nt z — 0.27 fm are due to ip[3i2] ~ Hi[z ~ 0) ~ 0. 
Calculations also show that although ^ the condition |s^| <C \s±\ is well fulfilled for 
most cases. 
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FIG. 2: The s± component of the spin density (in arbitrary units) in the ground state of ^°Mn, 
calculated for three selected cross sections at z=0.27fm (upper part) z=1.35 fm (middle part), and 
2=2.46 fm (lower part). Left panels show the contribution of the odd neutron right panels show 
polarization effect exerted by the odd neutron on the neutron core. For this case sj^ sj, see text 
for further detail. 
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a) 50Mn: 

7u[312]5/2«)v[312]5/2 
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b)5iMn: 

7u[312]5/2 



• 0.27fm 
o O.Slfm 
A 1.35fm 
V 1.90fm 
▼ 2.46fm 
n 3.03fm 

♦ 3.62fm 



c) 52Mn: 

7i[312]5/2(8)v[303]7/2 




40 80 120 160 40 80 120 160 40 80 120 160 

0(g^^(7t[312])^^(core)) 



FIG. 3: Product of |sj^^^^(r)||s^°^^(r)| reflecting the magnitude of the spin fields (in arbitrary 
units) versus the classical angle between these vectors 6, corresponding to their relative orientation, 
a) N=Z o-o nucleus 5°Mn, b) iV - Z = 1 odd- A nucleus ^^Mn, c) N - Z = 2 o-o nucleus ^^Mn. 
All points calculated at a fix value of the z coordinate and different values of the {x,y) coordinates 
are labeled by the same symbol, as indicated in the legend. 



Correlation between the spin field sj_ (r) due to the occupation of the [512]5/2 orbital 
by the valence proton and that of the core s^°''^(r) (polarization effect) is illustrated in 
Fig. El The figure shows the product |s^'^^^'(r)| ■ \s]^^'^{r)\ that reflects the magnitude of 
the spin fields versus the classical angle between these vectors 6 that reflects their relative 
orientation. The flgure clearly illustrates that in an o-o N=Z nucleus the core polarization 
is strongest and of almost purely ferromagnetic type. In an odd-A N—Z=±l nuclei the 
effect is quenched but remains to be of ferromagnetic type. Hence, the role of the spin flelds 
in these nuclei is maximal. In \N — Z\ > 1 nuclei the induced (core) spin fleld appears to 
be quenched even further and, additionally, not coherent with the valence-particle (s) spin 
field, particularly in o-o nuclei, as depicted in Fig. Efc. 

At least two very important conclusions can be drawn from this analysis, (z) The coher- 
ence of the spin fields in N^^Z nuclei may cause a strong polarization of the nucleus, (ii) 
The magnitude of the spin-field-induced effects is predicted to depend strongly on isospin. 
These two observations give an unique opportunity to resolve the strength of the spin fields, 
in particular, by using high spin states where spin fields are expected to be enhanced. 
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FIG. 4: Calculated energy differences for the terminating states AE't/^ = Eth[d~^^f^^^] — Eth[fy^] 
relative to the experimental data AE = AE'exp ~ ^Eth where AE(,xp are listed in Tab. d The left 
panel shows SHF calculations while right part illustrates calculations using l-ledf. 



B. The spin fields at the band termination 

Fig. shows the calculated energy differences for the terminating states AEth = 
['^3/2/7/2^] ~ ^th[fy2] relative to the experimental data AE^xp given in Tab. m i.e. the 
values of AE = AE^xp — AEth- The first striking observation stemming from this calculation 
is that all considered Skyrme forces systematically underestimate empirical data by at least 
10% . In case of the SkM* and SkP parameterizations the difference even exceeds 20-i-30%. 
The disagreement is unexpectedly large provided the structural simplicity of the terminating 
states. 

Let us further observe that the values of AE calculated using SLy4 and SLy5 forces, which 
are in general rather similar to AE obtained using SIII or SkO forces, increase rapidly in 
N—Z=l nuclei ^^V and '^^Ti. This result is related to the magnitude of the spin fields 
and the ferromagnetic-type polarization of the core which, as discussed in the preceding 
section and in Ref . |^ , are exceptionally large for Lyon- forces, see also Fig. [T] Since no 
enhancement of this type is observed in the data, this result clearly shows that an unified 
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description of the spin fields witliin tlie ledf tlieory is required. 

Hence, we state that the experimental data suggests a generalization of the s-ledf. 
However, our strategy is to introduce a minimal-type modification that pertains ultimately 
to the TO part of the S-LEDF only. More precisely, for SLy4, SHI, SkO, and SkM* we will 
change the spin fields i.e. the first two terms of the to part of the ledf (0) not affecting 
the local gauge invariance (jH)). For SkP, SkXc, and SLy5 on the other hand, we also slightly 
modify the to part of the tensor term. In this way we actually brake the local gauge 
invariance. However, our calculations show that this has a very small effect on the final 
results when compared to calculations using the Skyrme-force-induced Cj values. Thus, 
one can state that the local gauge invariance (jUI) is in fact preserved in our calculations. 

Our favorite unification scheme for the treatment of spin fields (called l-ledf) follows the 
one developed by Bender et al. js^. Let us recall, that in the l-ledf calculations we assume 
density independent coefficients defined through the Landau parameters 5^0=0.4, 5fQ=1.2, 
5'i=-0.19, and g[=Q.Q2 and set C^'^=Q. Such a simple treatment of the spin fields leads to 
a surprisingly consistent picture for various Skyrme-forces, irrespectively the difference in 
effective mass. Indeed, the results obtained for all forces, except of SkM* and SkP, essentially 
overly each other as shown in Fig.|3jD. Let us further observe that Ai? calculated with SkM* 
and SkP show almost a constant offset as compared to the other forces. We suspect that 
the effective mass or, equivalently, current independence of our results is directly related to 
the gauge invariance of the ledf. This point requires, however, further investigation. 

Since we do concentrate on N^Z nuclei, our calculations are essentially insensitive to 
changes in the isovector Landau parameters (7q and g[ in the wide range of their proposed 
values. Sensitivity of our predictions with respect to the isoscalar Landau parameter 
is illustrated in Fig. Fig. shows the average deviation from the data, AE, versus Qq. 
Note, that AE is minimal for ~ 0.4 -=-0.8, what is very close to the suggested value go=OA 
of Ref. 13]. Let us further observe that AE does not change sharply within the interval 
Ago = ±0.4 around the preferred value, but our analysis seem to rule out both negative and 
large positive go > 1.2 values of go. Moreover, it clearly shows that the ~10% discrepancy 
between the calculations and the data cannot be accounted for by further readjusting the 
Landau parameters, at least not within the analyzed unification scheme. 

Fig. shows the dependence of the standard deviation, ctas, which reflects the spread 
in AE, on go- Apparently the minimum is obtained for go ~ 0.8 -^ 1.2, i.e. well above the 
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FIG. 5: The average difference AE (AE = AE^xp — AEth) (upper part) and standard deviation a 
(lower part) between the data and the calculations versus Landau parameter qq. 



preferred value of gQ=QA. Let us observe, however, that almost all curves shown in Fig. 0] 
show clearly an increasing trend as a function of the reduced isospin = {N — Z) / A. Hence, 
part of the spread may merely reflect the isovector properties of the ledf, most likely, the 
isovector part of the £s-term. The relatively weak dependence of AE on Ta obtained for the 
SkO force seems to support this conclusion. Additional analysis strengthening this scenario 
will be given in the next section. 

V. THE SPIN-ORBIT TERM 

Within the SHF theory, the £s-potential takes the following form 

VLs{q.r) = -zW,(r)V x o", 
13 



(11) 



1 — I 1 r 



SkM* 



% o 



SkO 




SkXc 


• 






FIG. 6: The strength Wq (lower part) and the effective-mass scaled strength (upper part) of 
the isoscalar part of the spin-orbit SHF potential versus the effective mass. Filled dots denote 
parameterizations which do not include tensor densities while those including these terms are 
indicated by stars. See text for more detail. 



where 

W,{r) = ^WVpoir) + ^W'Vp.ir) - i [{t.x, + hx^)J{r) + {t^ - h)J^{r)] . (12) 

The vector spin density J{r) dependent terms appear to contribute rather weakly to the 
£s-potential. Hence, the magnitude of the £s-potential is determined essentially by the first 
two terms in Eq. p2|) . For the spherical limit the ^s-potential can be approximated by: 

V^L5(g,r)^ |l^ip^(r) + W^'ip;(r)|£s= |^p^(r)±^p;(r)|£s (13) 

14 



where pg = {pn + Pp) ' and p[ = {pn — Pp) ' are the radial derivatives of the local isoscalar and 
isovector densities, while Wq = W + and Wi = denote the isoscalar and isovector 
strengths, respectively. 

A direct comparison of the isoscalar strengths Wq of the £s-potential is given in Fig. IHt- 
Apparently, SLy4, Sly5, SkM*, and SIII forces have strong, while SkO, SkP, and SkXc have 
a weak £s-potential assuming, of course, that there are no drastic differences in the isoscalar 
density profile. The latter assumption should be rather well fulfilled in light nuclei, which 
are considered here. It is interesting to note that the conclusions stemming from a direct 
comparison of Wq are in complete contradiction to the results presented in Fig. |3] Indeed, 
according to our calculations Sly4, SLy5, SkO, SIII, and eventually also SkXc are expected 
to have similar £s-strength while it should be considerably stronger for SkM* and SkP. 

The question therefore arises of how to compare the strengths of the ^s-potential for 
different parameterizations. The problem appears to be related to non-local effects which are, 
within the SHF, absorbed into the kinetic energy term through the effective mass m*. The 
impact of non-localities on the £s-potential can be studied using the so called asymptotically 
equivalent wave function [3, 0| 

Hr) = \f^Mr). (14) 

This representation allows to rewrite the SHF equations in an alternative form with a bare 
mass m in the kinetic energy term, a state-dependent central potential Uq{ep_;r), and the 
effective- mass-scaled fs-potential (fT^ : 

nsM^'^['^pUr)±'^ri(r)]is. (15) 

The effective-mass-scaled isoscalar strengths Wq = are depicted in Fig. IHb- Note, 

that the classification of the £s-strength according to Wq agrees very nicely with our results 
shown in Fig. HI Indeed, the values of Wq are similar for Sly4, SLy5, SIII, and SkO and 
considerably larger for SkP and SkM*. 

This indicates that at least part of the observed discrepancy, AE, results from a too 
strong is-term. By reducing the strength with ~ 5%, the AE reduces by ~350 keV bringing 
it to an acceptable level of AE ~ 200 keV for most of the forces. In particular, for the case 
of the SkO force AE drops from 504 keV to 164 keV, while for Sly4 a 5% reduction of the 
£s-strength decreases AE from 560 keV to 189 keV. 
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Concerning the isovector £s-potential, the Skyrme forces which are discussed in the hter- 
ature can be divided into three major classes. The standard Skyrme force parameterizations 
assume that W = W' {Wi/Wo = 1/3), implying that Wg ~ Vr(2p^ + p'_g). The Sly4, Sly5, 
SkM*, SkO, and SIII are standard forces among those studied here. 

Non-standard Skyrme interactions with W ^ W were first studied by Reinhard and 
Flocard [38] in connection with isotope shifts in Pb nuclei. Consistency with experimental 
data led them to the parameterizations with W ~ —W or (H^i/ll^o = i-e- to an entirely 
different isovector dependence of the is-term Wg ~ W p'_^ as compared to the standard one. 
The study of the £s-term in neutron-rich nuclei by Reinhardt et al. seems further to 
corroborate this result. The so called SkO parametrization established in Ref. 3^ (and 
studied here) have even larger negative value of Wi with Wi/Wq ~ —1.3. 

The third type of the ^s-term which is considered in the literature in connection with the 
SHF approach, was introduced by Brown who uses W = Q (parametrization SkXc). In 
this case there is no isovector £s-term (W^i/VFo = 0) and Wg ~ W{p^ + Pp). 

As already discussed at the end of Sect. IIV Bl our calculations give certain preference for 
the SkO-induced l-ledf, since it minimizes the spread in AE, a^E- In particular, this 
result seems to speak in favor of a £s-potential with large negative isovector strength Wi. 
To reinforce this observation we have performed a set of calculations based on SkO-induced 
L-LEDF, but exploring different isovector dependence of the is-term, including the four 
possibilities discussed above Wi/Wq = —1.3,-1,0,1/3. In the calculations, the isoscalar 
strength Wq was kept constant and its value was reduced by 5% as compared to the original 
SkO strength. 

The calculated values of AE for these four variations of the SkO-induced l-ledf are 
shown in Fig. [3 Drastic change in the Wi/Wq ratio from the (near)original values -1.3,- 
1 to 0,1/3 clearly destroy the agreement to the data. Indeed, the spread, (Tab, increases 
from 113 keV and 136 keV {Wi/Wq = -1.3, -1) to 184keV and 180 keV {Wi/Wq = 1/3, 0), 
respectively, see insert in Fig.[7| Note, however, that since we deal with N^^Z, a fine tuning 
of the isovector terms cannot be achieved. 

Similar calculations, but with the SLy4-induced l-ledf additionally corroborate our 
conclusions. The change in Wi/Wq ratio from the (near) original values 1/3,0 to -1,-1.3 
improve the agreement to the data, as shown in the insert in Fig. Note also, that the 
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FIG. 7: The values of AE calculated using SkO-induced L-LEDF with four different parameteri- 
zations of the isovector £s-term, including the original Wi/Wq ~ — 1-3 strength and the modified 
strengths Wi/Wq = —1) 0, 1/3. In the calculations Wq was kept fix at 5% below its original value. 
The insert shows the dependence of the spread in AE, ctae, on Wi/Wq calculated for SkO-induced 
as well as SLy4-induced l-ledf. 

calculated spread, (Jae, is quantitatively very similar for both the SkO-induced and SLy4- 
induced l-ledf provided the isovector part of the is-term is similar. 
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VI. 2p-2h TERMINATING STATE IN ^^Sc AND THE ADDITIVITY OF Ip-lh EX- 
CITATIONS. 

The 2p-2h 0^^^/2/7/2^ terminating states can provide further insight into the to terms, 
the is potential, and into various polarization phenomena, since they allow to test self- 
consistent calculations against simple, intuitive additivity relations. Thus far such state is 
known, however, only in ^^Sc Although a single case does not allow for deep conclusions, 
the comparison between the data and the calculations, see Fig. |H1 clearly shows that some 
naively-understood additivity concepts do not always work. 

For example, the values of AE2p2h = ^-^2^2/1'' ~ ^-^2p2i calculated using Sly4, SLy5, 
SIII, and SkXc l-ledf are comparable or even smaller than the corresponding values of 
AE obtained for the terminating Ip-lh state in '^^Sc, see Fig. H] Consequently, calculations 
using a 5% reduced isoscalar £s-strength overestimate the data with exception of the SkXc 
L-LEDF for which one obtains an almost perfect agreement to the data. On the contrary, 
the additivity seems to work perfectly for the SkO-induced l-ledf. A reduction of the 
£s-strength by 5% gives, in this case, an excellent agreement to the data. 

The excitation energy of the 2p-2h terminating state in ^^Sc can be approximated by a 
sum of the energy of the two aligned 29/2+, [d'^/'2fj/2^]iiniax-i) states terminating at unfa- 
vored signature. Since there exist exactly two such states, their centroid energy (sum of 
their energies) should be fairly insensitive to the possible beyond-mean-field mixing effects. 
Hence, the additivity relation for the excitation energies should provide an independent 
characteristics of the ledf. 

Insert in Fig. |S1 shows the theoretically calculated difference AE^^^f^^ — AE2p2h where 
Ai?2p2/i denotes excitation energy of the terminating 2p-2h state while A£'2p2/!'*''^ denotes 
2p-2h excitation energy calculated by adding excitation energies of the two Ip-lh configu- 
rations given above. In general, the additivity works relatively well. There are however two 
exceptions: the SkP and, unexpectedly, SkXc l-ledf. The downward shift of SLy5 with 
respect to SLy4 as well as rather mediocre agreement for SkP and SkXc forces, may be is 
related to the tensor component that is present in these three forces and deserves further 
investigation. 

This short analysis indicates already the rich physics that can be addressed using 2p-2h 
terminating states. 
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FIG. 8: A£^2p2/i calculated for different Skyrme-force-induced L-LEDF. Open circles mark cal- 
culations using the standard ^s strength, while those using a 5% reduced ^s strength are la- 
beled by filled circles. The insert shows the theoretically calculated deviation from the additivity 
^E^p2h^^ - ^^tSh- The corresponding experimental value of AE^lf^^''^^ - A£;^p2? =0-252 MeV 
is marked by a solid line. 



VII. SUMMARY 

We have performed a systematic study of terminating states in the A~50 mass region 
using the self-consistent Skyrme-Hartree-Fock model and testing several parameterizations 

of the Skyrme force. The objective was to demonstrate that the terminating states, due 
to their intrinsic simplicity, offer an unique and so far unexplored opportunity to study 
different aspects of effective NN interaction or nuclear local energy density functional within 
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FIG. 9: The average deviation IS.E for N^Z (open dots) and for N=Z (filled dots) nuclei. The 
latter case includes data in ^^Ti 19] where ^^[/^/g] 12+ =8-038 MeV, E[d~l^f:i'l:^\^- =U 2m MaY , 
and A£;exp=5.331 MeV, and ^'^V :39'] where E[fJ}j^]^^+=^AMMQ\, E[d-lj:}^^\'j - =l?,.m<d MeY , 
and A£'e2;p=5.145 MeV. All calculations were done using l-ledf. 



the self-consistent approaches. 

It is shown, that the Skyrme-force parameterizations used in our work, including Sly4, 
Sly5, SkO, SIII, SkXc, SkP, and SkM*, give rise to a rather mediocre (for SkP and SkM* 
even unacceptable) agreement with the data even for such a seemingly simple observable 
like the energy difference between the aligned fy^ and c^j^^j/T/t^ states. 

It is further demonstrated that a simple unification of the spin fields according to the 
scheme proposed in Ref. 3| leads to an unified description of the data for Sly4, Sly5, 
SkO, SIII, SkXc, i.e. for very different parameterizations. This result seems to indicate 
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the importance of the local gauge (and Galilean) invariance of the local energy density 
functional. The remaining discrepancy of ~500 keV (see Fig. ^ which is still of the order of 
~10% cannot be reduced by further readjusting the spin fields. 

It is also shown, that the observed disagreement between theory and experiment for differ- 
ent parameterizations correlates nicely with the values of the effective mass scaled isoscalar 
strength of the is-term for these parameterizations. Hence, a part of this discrepancy can, 
most likely, be ascribed to a too strong isoscalar is-term. Reduction of the isoscalar is- 
strength by 5% reduces the discrepancy well below the 5% level. Moreover, our calculations 
suggest that the spread in AE can be further reduced by adopting the non-standard param- 
eterizations of the is-term with a strong negative isovector strength Wi/Wq < — 1. 

Finally, let us point out that there is a large difference in AE calculated in N^Z and 
N=Z nuclei, see Fig. El For Ny^Z our l-ledf approach systematically underestimates the 
data while the opposite is true for N=Z nuclei. The offset between the two curves is of 
the order of ~1 MeV. This result leaves us with an extremely important question: Does 
this offset indicate a breakdown of the standard mean-field in iV~Z nuclei and a need for 
substantial configuration mixing even at the terminating states j^l and what is the possible 
source of such mixing? 
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